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Abstract 

We study the complete time evolution of scalar fields propagating in space-times of higher di- 
mensional Lifshitz Black Holes with dynamical critical exponent z = 2, obtained from a theory 
including the most general quadratic curvature corrections to Einstein-Hilbert gravity in D di- 



5-H 

mensions. We also computed the quasinormal spectrum after performing a numerical integration 
and solving exactly the Klein-Gordon equation obeyed by the massive scalar field. We found that 
quasinormal modes are purely imaginary for all dimensions. 
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I. INTRODUCTION 



The study of black hole perturbations has been of interest since the original work of 
Regge and Wheeler [l]. The consideration of the fluctuations of the background space-time 
geometry is related to experiments for the detection of gravitational waves from supernovae 
or from coalescence of binary neutron stars, that are supposed to be candidates for the 
formation of black holes. During this process, gravitational waves are emited, a process that 
at late times is dominated by the fundamental quasinormal mode 

Apart from gravitational fluctuations, it is useful to consider perturbations of matter 
fields with different spin weights outside the Black Hole space-times. It is know that the 
time evolution of such perturbations can in general be divided in three stages: a short 
period of initial outburst of the perturbation, followed by a long period characterized by 
proper damping quasinormal oscillations, and at late times, depending of the asymptotic 
of the space-time at spatial infinity, a stage in which quasinormal modes are suppressed by 
power-law or exponential tails js-7|. 

Quasinormal modes are entirely fixed by the geometry of the background Black Hole 
space-time and are independent of the initial conditions on the perturbations. There are 
similarities with the usual normal modes of closed systems, but the open nature of the physi- 
cal region related to the fact that the fluctuations can be absorbed by black hole or radiated 
to infinity, lead to frequencies that in general are complex, and do not form a complete 
set, dominating the evolution over a long but well defined period of time. As quasinormal 
frequencies depends on the parameters of the black hole geometry, their computations can 
serve to estimate such parameters js]. 

On the other hand, the well known relationship between Anti de Sitter (AdS) bulk space 
gravitational problem in D dimensions and Conformal Field Theory (CFT) at the (D — 1)- 
dimensional border, also known as AdS / CFT relation, extrapolated its original applications 
to be studied in many areas, including phase transitions in the CFT side and generalizations 
of gravity in the bulk. Quasinormal modes are also important in this context 9KLll]. Related 
to this idea, a black hole in asymptotically AdS space-time is dual to a strongly interacting 
system defined at the AdS boundary, and the imaginary part of the fundamental quasinormal 



mode is proportional to the thermalization time in the boundary theory 



12 



13(]. The 



gauge/gravity duality has been applied to study a great variety of physical systems, from 
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quark-gluon plasmas to holographic superconductors 

Recently, gravity solutions with strong asymmetry of time with respect to space coordi- 
nates have been obtained, both in the severely modified Horava Lifshitz gravity [16] as well 
as in Einstein gravity with higher order curvature corrections 17]. 

Several theories present anisotropic scaling when one compares time with space coordi- 
nates. Such rescaling appears frequently when higher powers of curvature are present in the 
gravity action. The so called D-dimensional Lifshitz backgrounds have the form 

ds 2 = -^dt 2 + —dr 2 + dx a dx a , (1) 

Jj2z rp2 J^2 / ^ ' v ' 

i=l 

where x a are the components of a (D — 2)-dimensional coordinate vector. The geometry 
described by the above line element admits the anisotropic scaling symmetry 

t H> X z t , r ^ A~V , x a H> Xx a , (2) 

which is precisely part of the isometry group present in their non-relativistic duals, that are 
considered on the (D — l)-dimensional boundary at infinity. 

This paper is the first of a series in which a complete study of the time evolution of per- 
turbations around Black Hole solutions whose asymptotic is the space-time (CE]) is performed. 
We will consider scalar perturbations, that are important not only due to its simplicity, but 
for the implementation of gauge/gravity duality ideas. Our motivation is basically concerned 
with the extension of the AdS / CFT correspondence to other areas of physics, including fur- 
ther condensed matter models. Within the context of non-relativistic physics, gauge/gravity 



ideas were recently addressed, as in references 



18] and 19] , where the authors proposed grav- 



ity duals to non-relativistic systems, and in [20], where new gravitational duals to Lifshitz 
fixed points with scale invariance but without Galilean invariance were presented. 

Previously, the three dimensional case had been considered, the so called New Massive 



Gravity 21], and it has been conjectured that such a theory is dual to the KdV equation in 



22] . It will be encouraging if one is able to find similar relations 



two dimensional space-time 
in higher dimensions. 

The paper is organized as follows: Section II briefly presents the line element of higher 
dimensional asymptotically Lifshitz black holes obtained considering the most general 
quadratic corrections to Einstein gravity. In Section III we obtain a general wave equa- 
tion suitable to analyze the propagation of scalar fields in the background geometry of the 



higher dimensional black holes, whereas Section IV is devoted to the numerical solution of 
the evolution equations, and the computation of the complex quasinormal frequencies us- 
ing two methods: fitting time domain data and the Horowitz-Hubbeny approach. We also 
present, in Section V, the exact results for quasinormal frequencies after solving the equation 
of motion for the scalar field in the general .D- dimensional Lifshitz backgrounds. Finally, we 
present our conclusions. 



II. D > 5 ASYMPTOTICALLY LIFSHITZ BLACK HOLE FAMILIES FROM 
QUADRATIC CORRECTED GRAVITY 

In this section we briefly present three different classes of asymptotically Lifshitz black 
hole solutions which correspond to different ranges of the dynamical exponent z • As was 
shown in [l^] below, there exist at least one family of black hole solutions for any generic 
value of the dynamical exponent z. 

In the following we will consider the action: 

S = J dPxy/=g (R-2A + ^R 2 + &R, V R^ + &R, var R^ T ) . (3) 

that includes the most general quadratic-curvature corrections to D-dimensional gravity, 
and we will focus only in the higher dimensional cases with D > 5. 
The field equations that arise from the action are jl7| : 



G>, + Ag^ + + Ag^UR + i^V^i? + $ A R^ VT R ar + 



■fJ,V 



+ 2^3 {R^arRv^ — 2R tM7 R") — - [ilR 2 + £,2RarR aT + ^RuT-fdR^ 15 ) g^u = . (4) 

where the coefficients k — 1, ...,4 in front of some geometric factors are functions of 
those appearing accompanying the quadratic terms in the action (j3j) and are given by: 

tfi = 6 + 4£ 3 (5) 

02 = ^(4£i + 6) (6) 

tf 3 = - (26 + 6 + 2£ 3 ) (7) 

tf 4 = 2 (6 + 26) (8) 
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The above field equations have solutions whose line elements are particular cases of those 
describing a D-dimensional manifold M, and that can be written locally as a warped product 
of the form 

M D = 2 xlC D - 2 , (9) 

where we assume that M m is a lorentzian and KP~ 2 an Einstein manifold. Introducing 
coordinates z^ = (y a ,x l ) in A4 D , where y a and x l are coordinates defined in the manifolds 
O 2 and JC D ~ 2 , respectively, the full metric describing the geometry of the manifold Ai D can 
be written as 

ds 2 = Q^dz»dz» = gab (y a )dy a dy b + w 2 (y a )da 2 D _ 2 , (10) 
where w(y a ) is an arbitrary function defined in O 2 and 

da 2 = ^ij{x)dx l dx^ (11) 

is the metric on KP~ 2 . 

A common realization of (fTUj) is that in which we chose radial and time coordinates 
(y 1 , y 2 ) = (r, t) in O 2 and describe the geometry of M. D by the line element 

ds 2 = -A(r)dt 2 + B{r)dr 2 + C(r)dE^_ 2 (^) . (12) 

As an example we can mention the Lifshitz spacetimes (jTJ), that are solutions of (JH) for a 
general value of the dynamical exponent z in any dimension, in such a way that they exist 
for fixed values of the cosmological constant and the parameter given by 

1 / 2 , w s i(D-3)(D-A)z(z + D-2)B 3 \ , . 
A = -^[2z 2 + (D-2)(2z + D-l)-^ ^ ,(13) 

L 2 - 2 [2z 2 + (D - 2)(2z + g - 1)] & - 4 [z 2 — (D — 2)z + 1] /3 3 
^ 2 " 2^ + D-2) ' (M) 



generalizing the results for the D = 3 New Massive Gravity obtained in Ref. 26| for /3 3 = 



and /?2 — — (8/3)/5i = — 1/m 2 . The formula ([T]) can be written in the form (|12p with 
A(r) = gl, S(r) = C(r) = g and 7i ,(^) = diay[l, 1, 1], i, j = 1, ..d - 2. 

If the line element ( fTUj) describes a black- hole space-time (that is not the case of (fH)), 
then /C n describes the structure of a spatial section of its event horizon. If K, n is a constant 
curvature space with sectional curvature K, then = represents a flat space, if = 1 a 
spherical manifold and K = — 1 a hyperbolic one j^ . 



For D > 5 there are three different Lifshitz black hole families. The geometry of all these 
families are described by (fl2]l with C(r) = and dT? D _ 2 (x l ) = dx a dx a , a = 1, ..D — 2 
represents a flat euclidean (D — 2)-dimensional manifold with K — 0. The first family of 
solutions describes asymptotically Lifshitz black holes for the dynamical exponent z > 2 — D, 
and its line element is given by ( 1121) with 



. . r 2z ( ML^ D ^I 2 \ , . L 2 ( ML^+n-^Y 1 , . 

^) = Jir( v 1 - r ( 8+g -2)/2 j ' B ( r ) = -2{ 1 - r (z + D-2),2 ) ( 15 ) 

The dynamical exponent z parameterizes the coupling constants of the theory. This solution 
is defined only for D > 5, and there exists a conformal limit z = 1 in the form of an 
asymptotically AdS black hole that coincides with a solution of the Einstein- Gauss-Bonnet 
gravity with a fine-tuned coupling constant. 

The second family of asymptotically Lifshitz black holes are obtained for z > 1, and are 
of the form (Tl2l with 

. . r 2z ( ML 2 ^\ /N L 2 ( ML^-^Y 1 , . 

and is compatible with the coupling constants given by 
(z-l)[z 2 -Dz-(D -!)(£> -2)] 



A 



2L 2 (z- J D) 



/3i = L 2 3(D - 1)(D - 2)^ 3 - (2D 3 - 2D 2 - UD + 20)z 2 + (3D 3 - IAD 2 + 19D + 10)^ 

+ ( J D + 2)( J D-4) / 2(D-2)(D-3)(D-4)«(2;- 1)(« - D)(3« + D - 4) 

[6(D - 2)z 2 - (D 2 -3D + 8)z - 2(D - 4)] 



02 = —L(D — l)(2z — D — 2) 



[2(D - 2)(D - 3)(D - A)z(z - l)(z - D)(3z + D - 4)] 

L 2 (D -l)(2z- D -2) 



133 ' A(D - 3)(D - A)z(z - D) ' 
In similarity with the above case for the z > 2 — D black hole family, the solution exists 
only for higher dimensions D > 5 and a conformal limit z — 1 is absent. 

In reference [l?| the authors also found a third family of Lifshitz black holes caracterized 
by a negative dynamical critical exponent z = —\z\, with metric given by ffT2l) with 



with the corresponding coupling constants parameterized as usual by z. Again this family 
of black holes is defined only in higher dimensions D > 5, and due to the presence of a 
negative dynamical critical exponents, it has no conformal limit z — 1. 

In the following we will be interest in the second family of higher dimensional Lifshitz 
black holes discussed above and specifically those solutions that are caracterized by the 
dynamical exponent z = 2. 



III. CAUSAL STRUCTURE OF THE z = 2 LIFSHITZ BLACK HOLE 



In this section we study the causal structure of z = 2 asymptotically Lifshitz black holes. 
Considering for simplicity the five dimensional case and taking L = 1, we can write the line 
element as 

ds 2 = -r 4 ( 1 - -t \ dt 2 + — -r^- — ^dr 2 + r 2 (dx 2 + dx\ + dxV) . (18) 
Regarding to the causal structure, the above metric is very similar to the z = 3 NMG black 



hole [2l|. There is a single regular event horizon at r + = \/M, where M is an integration 
constant related to de ADM mass. Also, the metric exhibits a spacetime singularity at r = 0, 
as we can see from the behavior of the Kretschmann invariant 

R,u. P R^ ap = ~ (3r 4 - 9r 2 + r 2 + 17r 4 ) . (19) 

As r — > we have 

We also notice that as r — > r + 

RiivapR^^ 9 — > —88 . 

Therefore, from these facts we can say that there is a genuine physical singularity at r = 
and a regular surface at r = r + . Further details can be extracted from the Penrose- Carter 
Diagram. In order to construct such a diagram, we have to remove the singularity in the 
coordinate patch at r = r + through the null coordinates 

= e r»(t+r.) ? V = _ e -r£(t-r.) ? 
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where r* is the tortoise coordinate given by 

r*(r) = In ^1 

The metric f lT8|) can be rewritten as 

ds 2 = - j d£W + r 2 (cfe 2 + + (fa 2 ;) , 

which is obviously regular at r = r+. 

Performing another change in the coordinates, that is, 

T = arctan(CZ) + arctan(V), X = arctan(?7) — arctan(V) , 

we construct the Penrose-Carter diagram (1T1). 



r=oo 




r=oo 



Figure 1: Penrose-Carter Diagram for five dimensional z = 2 Lifshitz Black Hole. The dashed lines 
refer to the null physical singularity at r = covered by a regular event horizon at r = r+. 



The diagrams show us the null nature of the singularity and the AdSdike spatial infinity 
at r = oo. From this, we can see that the causal structure of z = 2 Lifshitz black holes are 
essentially the same as the ordinary AdS black holes. Thus, one can apply the gauge/gravity 
duality ideas in the Lifshitz case, in a pattern similar to the AdS case. 
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IV. PERTURBATION EQUATIONS 



In the space-time (fTUj) the Klein- Gordon equation for a minimally coupled massive scalar 
field (biz^) with mass m can be written as 



t E ^{z v ) = (□ - m 2 )ct){z v ) = . 



(20) 



We can write the Laplace-Beltrami operator K E in (Jz0|) as □ = V V M = □ + ^D a w ■ D a + 
w~ 2 A, where □ = D a D a and A = D l Z)j are the Laplace-Beltrami operators in Af m and /C n . 
In the above formulas D a and -D* are the usual covariant derivatives in M m and /C n . 
Introducing scalar harmonics §fc s (a; 1 ) in /C n , that are eigenfunctions of A defined by 



(A + k 2 E )S kE = 0, 
we can separate variables to obtain 

Putting fl22|) in fl20l) . we obtain the equation for (f>(y a ,x l ) on A/"'' 

□y, kB (^) + WgW(v°) = mW?/ a ) 
where = W^ E \y a ) is given by the expression 



(21) 



(22) 



(23) 



k B 



— ^ — -(Dw) z - -wDw - k E 



(24) 



Above, are the eigenvalues of the scalar harmonics. Equation ( 12 3 p describes the propaga- 
tion of scalar perturbations over a generic (m + n)-dimensional space-time Ai. Considering 
now the particular case of a D-dimensional space-time given by (fl2"|) . with dY, D ._ 2 (x' 1 ) repre- 
senting a flat euclidean (D — 2)-dimensional manifold, then (123]) results in 



with = (t, r) and 



i r(^-2)(4-^) F(r) _^2 G(r) _, 2 



(25) 



C(r) 



4 



G(r) = V ^)D (2) ( V ^)) 



(26) 
(27) 
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where CU) denotes the Laplace-Beltrami in the two-dimensional orbit space with coordinates 
(y 1 ,?/ 2 ) = (£,r). Using the explicit metric components in the orbit space in (}25l)- (l2T{) we 
obtain 



2,„(^) 



A(r) \ gyjfj 1 / A 
B(r)J dr 2 + 2 If? 



d 
dr 



In 



(E) 



dr 



+ A 



W^\r)-m 2 



(28) 



The functions F(r) and G(r) in (r) are given by 

2 



1 J (D-2)(6-D) fdCy D-2 



C{r) 



16BC 



dr J 



8B 



-In | i 



dr 



B 



dC 



+ 



1 d 2 C 



dr 2B dr 2 



(29) 



r- It maps the 



We introduce in f l28l) the tortoise coordinate r* defined by , , ,, ,.. 
physical region r G (r + ,oo) into r* G (— oo,0). Dropping the indices k and (E) we finally 
obtain for the evolution equation 

d 2 (p(t,r) dip(t,r 



8t 2 



where 
V(r) 



B 



16 V A 

A 



AH{ - 8Hi 



dr \B 



+ 



A 



dr 2 



+ 3 



V(r)<p(t,r) 



(30) 



d (A 



dr\B 



4 



A d 2 (A 



Bdr 2 \B 



A d TT 



W^\r) -m 2 



(31) 



plays the role of effective potential. 



V. TIME EVOLUTION OF PERTURBATIONS 



Equation fl30|) can be solved numerically by finite difference method. Taking t = t + kAt 
and r* = r* + jAr*, substituting in (|30|) and rearranging the terms we obtain the difference 
equation 

A+2 / A f 2 \ 

ti+i = + ^1 K +1 + ^i" 1 ) + (2 - 2 ^ - Vj j ft (32) 

As initial conditions we choose the static gaussian 

V(r„,to) = Ae-^*" 6 ) 2 , J^(r.,*) |t=t = (33) 

We also impose Dirichlet conditions at the AdS boundary ip(r* = 0,t) = 0. The Von 
Neumann stability condition applied to the above difference equations results in the relation 

At 2 . o /a\ At 2 



Ar 2 



sin 2 (f ) + ^-V(a) < 1 
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Figure 2: Semi-logarithmic plots of the time-domain evolution of m = and m = 
bations in Lifshitz black holes for various dimensions. The parameters are M = 
from top to bottom curves D = 5, 10. 



1 scalar pertur- 
L = k = 1 and 



Then, defining V m3iX as the largest value of Vj in the numerical grid, our numerical solution 
will be stable if 



At 2 At 2 



Ar 2 + — V — < 1 ( 34 ) 
Now taking At = \Ar if and choosing a sufficiently small time step such that V max < 
we assure the stability of our numerical codes. In all the calculations performed by us we 
verified the fulfillment of this stability condition. 

As Figures[2]l5]shows, the time evolution of scalar fluctuations in z = 2 Lifshitz black holes 
can be divided in two stages. The first one is a transient response of the field due to its initial 
interaction with the black hole, and is strongly dependent on the initial conditions for the 
perturbation. This transient stage is followed by purely damped quasinormal modes, without 
oscillations, in contrast with the usual behavior of asymptotically AdS black holes, where 
quasinormal oscillations caracterized by frequencies with well defined real and imaginary 
parts dominate all the evolution after the transient phase. 

The above picture is valid for all dimensions, and values of the parameters M, L, m and 
k that characterized the background space-time and the scalar field fluctuations. As the 
dimensionality of the space-times increases, the amplitude of massless and massive scalar 
fluctuations decays more quickly. Then, Lifshitz black holes becomes more rigid for higher 
dimensions. The same behavior is observed for scalar fluctuations with increasing wave- 
number k, for all dimensions. For fixed dimensionality of the background, s-wave fluctuations 
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Figure 3: Semi-logarithmic plots of the time-domain evolution of scalar perturbations in Lifshitz 
black holes for various dimensions, for M = L = m = 1 and from top to bottom curves, k = 0, ...6. 

are more longer lived. For fixed k the higher the dimension of the space-time, the lower the 
decay rate of the perturbation. 

When the mass of black hole M and that of the scalar field fluctuation m change, we 
observe, for D and k fixed, an increase of the decay rate of the perturbation with the 
increase of those parameters. Then, the propagation of scalar perturbations outside higher 
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5 10 15 20 5 10 15 20 

t t 

Figure 4: Semi-logarithmic plots of the time-domain evolution of scalar perturbations in Lifshitz 
black holes for various dimensions, for k = L = m = 1 and from top to bottom curves, M = 1, ...6. 

dimensional Lifshitz black hole exhibits a very simple behavior. 

All of the above results are confirmed after the calculation of the quasinormal frequencies 
of the scalar perturbations in the systems considered. The quasinormal modes are solutions 
of equations ( 1361) or ( 1301) with the specific boundary conditions requiring purely outgoing 
waves at spatial infinity and purely incoming waves on the event horizon r + = a/ML. 

In practical numerical calculations, the time evolution data obtained solving the evolution 
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Figure 5: Semi-logarithmic plots of the time-domain evolution of scalar perturbations in Lifshitz 
black holes for various dimensions, for M = L = k = 1 and from top to bottom curves, m = 0, ...6. 

equations mix the fundamental quasinormal mode and the higher overtones. As higher 
overtones decay more quickly that the fundamental mode, a common recipe is to fit the 
numerical data far away from the beginning of the quasinormal ringing, to be sure that 
higher overtones are suppressed. 

Another numerical method to determine the quasinormal frequencies for AdS space-times, 
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that can be easily adapted to asymptotically Lifshitz space-times is the Horowitz-Hubbeny 



method 



12|. 



Tables I to VI show the numerical results obtained for the quasinormal frequencies of 
higher dimensional Lifshitz black holes with dynamical exponent z — 2, using the above 
mentioned numerical methods. We also include the exact results that can be obtained 
solving analytically the evolution equation, as we will show in the next section. 

As we can observe from the tables, there exist a very good coincidence between the 
analytical results and the numerical ones. The Horowitz-Hubbeny method appears to be 
the best to perform numerical calculations and the time domain integrations lead to little 
worse results, due to the presence of numerical errors and precision problems in a practical 
calculation. However, if an analytical formula is absent, the Horowitz-Hubbeny method is 
very difficult to implement, because we need an approximate knowledge of the quasinormal 
frequency to complete the calculation, a fact that indicates the importance of develop time 
domain integration. 

It is quite surprising that the evolution of perturbations is extremely simple. Indeed, the 
decay is basically that of an exponentially decaying tail, which means no real part of the 
frequency. This sounds like no quasi normal decay. This may be conjectured as a normal 
phase, or no superconductor at the border. In figure ([2]) one sees that the perturbations 
boringly die out leaving no sound of the Black Hole, just showing that perturbations spread 
and disappear. This result is independent of the dimensionality of space time and spin 
number. Therefore, the structure of Lifshitz Black Hole does not seem to lead to any 
interesting critical phenomena at the border. 

Moreover, the incredible similarity of decays, almost independently of dimension or m- 
value, with strikingly small differences in the values of the constants, shows that the z = 2 
Black Holes are much more stable than the Einstein gravity counterparts. This leads again 
to the conjecture that the CFT dual never undergoes phase transitions. Actually, this result 
goes in the same line as the three-dimensional New Massive Gravity, whose CFT dual is 
conjectured to be an integrable model, namely, KdV equation [22 1. 
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M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.677057* 


-1.677033* 


-1.677033* 


2 


1 


1 


-3.218724* 


-3.218659* 


-3.218659* 


3 


1 


1 


-4.760508* 


-4.760286* 


-4.760286* 


4 


1 


1 


-6.302442* 


-6.301912* 


-6.301912* 


5 


1 


1 


-7.844577* 


-7.843538* 


-7.843538* 


1 





1 


-1.541628* 


-1.541626* 


-1.541626* 


1 


2 


1 


-2.083253* 


-2.083253* 


-2.083253* 


1 


3 


1 


-2.760615* 


-2.760286* 


-2.760286* 


1 


4 


1 


-3.709838* 


-3.708132* 


-3.708132* 


1 


5 


1 


-4.917540* 


-4.926791* 


-4.926791* 


1 


1 





-1.571455* 


-1.571429* 


-1.571429* 


1 


1 


2 


-1.952045* 


-1.952027* 


-1.952027* 


1 


1 


3 


-2.325160* 


-2.325145* 


-2.325145* 


1 


1 


4 


-2.749181* 


-2.749172* 


-2.749172* 


1 


1 


5 


-3.200247* 


-3.200247* 


-3.200247* 



Table I: Scalar quasinormal frequencies in 5 -dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 



VI. EXACT QUASINORMAL FREQUENCIES 

Equation (1251) can be analytically solved, offering us the possibility of find exactly the 
quasinormal frequencies of Lifzhitz black holes considered in this report. 
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M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.720780* 


-1.720759* 


-1.720759* 


2 


1 


1 


-3.321490* 


-3.321392* 


-3.321392* 


3 


1 


1 


-4.922364* 


-4.922025* 


-4.922025* 


4 


1 


1 


-6.523469* 


-6.522657* 


-6.522657* 


5 


1 


1 


-8.124882* 


-8.123290* 


-8.123290* 


1 





1 


-1.600632* 


-1.600633* 


-1.600633* 


1 


2 


1 


-2.081199* 


-2.081139* 


-2.081139* 


1 


3 


1 


-2.682273* 


-2.681771* 


-2.681771* 


1 


4 


1 


-3.523343* 


-3.522657* 


-3.522657* 


1 


5 


1 


-4.565130* 


-4.603787* 


-4.603787* 


1 


1 





-1.571455* 


-1.571455* 


-1.571455* 


1 


1 


2 


-1.977047* 


-1.977082* 


-1.977082* 


1 


1 


3 


-2.335233* 


-2.335205* 


-2.335205* 


1 


1 


4 


-2.750061* 


-2.750000* 


-2.750000* 


1 


1 


5 


-3.659501* 


-3.195887* 


-3.195887* 



Table II: Scalar quasinormal frequencies in 6- dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 

In order to compute the quasinormal frequencies that dominated at intermediate times, 
we assume for the function (p(t, r) in equation ( J28l) the time dependence 

ip(t,r) = Z{r) exp(— iut) (35) 

Then, for z = 2 higher dimensional Lifshitz black holes whose metric tensor have the com- 
ponents given by ( |T751) . the scalar field amplitude Z(r) satisfies the equation 

r 2 (r 2 - ML 2 ) Z"+r \(D + 1) r 2 - (D - 1) ML 2 ] Z'- ( L 4 k 2 + L 2 r 2 m 2 - ^-^ ) Z = 

\ r 2 — ML 2 J 

(36) 

were we denote the radial derivatives of a function f(r) as /' and drop from (p in (|28|) the 
labels k and (E). 
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M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.754315? 


-1.754301? 


-1.754301? 


2 


1 


1 


-3.400939i 


-3.400845?' 


-3.400845i 


3 


1 


1 


-5.047710i 


-5.047389i 


-5.047389i 


4 


1 


1 


-6.694699i 


-6.693934? 


-6.693934? 


5 


1 


1 


-8.341979?' 


-8.340478i 


-8.340478i 


1 





1 


-1.646544? 


-1.646544? 


-1.646544? 


1 


2 


1 


-2.077628i 


-2.077573i 


-2.077573i 


1 


3 


1 


-2.616515? 


-2.616360i 


-2.616360i 


1 


4 


1 


-3.370879?' 


-3.370661? 


-3.370661? 


1 


5 


1 


-4.325313i 


-4.340478i 


-4.340478i 


1 


1 





-1.666672? 


-1.666667?' 


-1.666667?' 


1 


1 


2 


-1.993564? 


-1.993813i 


-1.993813i 


1 


1 


3 


-2.336972?' 


-2.336953i 


-2.336953i 


1 


1 


4 


-2.741867i 


-2.741864i 


-2.741864i 


1 


1 


5 


-3.182104? 


-3.182091? 


-3.182091? 



Table III: Scalar quasinormal frequencies in 7 -dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 



Introducing in the above equation the new variable ((r) defined as 

ML 2 



C = i 



such that C( r +) — and C(°°) — 1) we obtain 
C(C-l)Z"(C) + \[(D-A)C + 2] Z'(C) + L 2 



*„2 ^ 2 (C-1) rn 2 M 2 

+ — ^ — - - 

Proposing a solution in the form Z(() = ( a (l — ()PF((), where 

iLui 



(37) 



Z(C) = (38) 



a = 



2M 



(5 = \D + \^Am 2 L 2 + D 2 
4 4 



(39) 
(40) 



18 



M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.7807875* 


-1.780776? 


-1.780776* 


2 


1 


1 


-3.464048*' 


-3.463956*' 


-3.463956* 


3 


1 


1 


-5.147450*' 


-5.147135* 


-5.147135* 


4 


1 


1 


-6.831065*' 


-6.830314* 


-6.830314* 


5 


1 


1 


-8.514963*' 


-8.513494* 


-8.513494* 


1 





1 


-1.683179? 


-1.683179* 


-1.683179* 


1 


2 


1 


-2.073610*' 


-2.073567*' 


-2.073567*' 


1 


3 


1 


-2.561607*' 


-2.561553*' 


-2.561553*' 


1 


4 


1 


-3.244799*' 


-3.244732*' 


-3.244732* 


1 


5 


1 


-4.111895* 


-4.123106* 


-4.123106* 


1 


1 





-1.700032* 


-1.700000* 


-1.700000* 


1 


1 


2 


-2.005357*' 


-2.005092* 


-2.005092* 


1 


1 


3 


-2.333347*' 


-2.333333*' 


-2.333333*' 


1 


1 


4 


-2.727540*' 


-2.727543* 


-2.727543* 


1 


1 


5 


-3.161256* 


-3.161250* 


-3.161250*' 



Table IV: Scalar quasinormal frequencies in 8 -dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 

we can transform fl38l) into an hypergeometric equation 

C(l - + [c - (1 + a + b)Q F'(C) - abF(() = , (41) 

whose general solution can be written as 

Z (C) = C (1 - Cf [Ax 2 F t (a, 6, c, () + A 2 C^ c 2 *i (6 - c + 1, a - c + 1, 2 - c, 0] , (42) 
where 2-F1 (C) is the hypergeometric function, and Ai and A2 are integration constants. The 
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M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.802098* 


-1.802177* 


-1.802177* 


2 


1 


1 


-3.515213i 


-3.515213* 


-3.515213* 


3 


1 


1 


-5.228248* 


-5. 228254* 


-5.228254* 


4 


1 


1 


-6.941137* 


-6.941296* 


-6.941296* 


5 


1 


1 


-8.654337* 


-8.654337* 


-8.654337* 


1 





1 


-1.713041* 


-1.713041* 


-1.713041* 


1 


2 


1 


-2.072675* 


-2.069562* 


-2.069562* 


1 


3 


1 


-2.551524* 


-2.515213* 


-2.515213* 


1 


4 


1 


-3.139086* 


-3.139124* 


-3.139124* 


1 


5 


1 


-3.941332* 


-3.941296* 


-3.941296* 


1 


1 





-1.727273* 


-1.727273* 


-1.727273* 


1 


1 


2 


-2.012838* 


-2.012756* 


-2.012756* 


1 


1 


3 


-2.012838* 


-2.326399* 


-2.326399* 


1 


1 


4 


-2.709162* 


-2.709208* 


-2.709208* 


1 


1 


5 


-3.135388* 


-3.135421* 


-3.135421* 



Table V: Scalar quasinormal frequencies in 9- dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 

arguments of the hypergeometric function are 

a = a + p - - + j^^J{ D ~ 2 ) 2 M 2 - AL 2 (Mk 2 + u 2 ) , (43) 
b = a + f3 - - - JL yJ(D - 2) 2 M 2 - AL 2 (Mk 2 + u 2 ) , (44) 
c = 1 + 2a . (45) 

Taking into account (14"5]) we can write (142]) as 

Z(C) = AiC a (l-Cf 2 i 7 \(a,6,c,C) + A 2 C Q (l-C) /3 2j Pi(6-c+l,a-c+l,2- C) C) , 

(46) 

Ingoing boundary condition at the horizon ( = of the Lifshitz black holes fixes A2 = 0, 
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M 


k 


m 


Time integration 


Horowitz-Hubbeny 


Analytical 


1 


1 


1 


-1.819814* 


-1.819804* 


-1.819804* 


2 


1 


1 


-3.557720*' 


-3.557627*' 


-3.557627i 


3 


1 


1 


-5.295766*' 


-5.295451* 


-5.295451* 


4 


1 


1 


-7.034023* 


-7.033274*' 


-7.033274*' 


5 


1 


1 


-8.772563*' 


-8.771098*' 


-8.771098*' 


1 





1 


-1.737823* 


-1.737823* 


-1.737823* 


1 


2 


1 


-2.065770*' 


-2.065745i 


-2.065745i 


1 


3 


1 


-2.475657i 


-2.475647i 


-2.475647i 


1 


4 


1 


-3.049581* 


-3.049510*' 


-3.049510*' 


1 


5 


1 


-3.788032*' 


-3.787333*' 


-3.787333*' 


1 


1 





-1.750000* 


-1.750000* 


-1.750000* 


1 


1 


2 


-2.018032*' 


-2.017985*' 


-2.017985*' 


1 


1 


3 


-2.317541* 


-2.317532*' 


-2.317532*' 


1 


1 


4 


-2.688471* 


-2.688476*' 


-2.688476*' 


1 


1 


5 


-3.106290*' 


-3.106289*' 


-3.106289*' 



Table VI: Scalar quasinormal frequencies in 10 -dimensional L = 1 Lifshitz black holes for various 
values of M , m and k. Results of calculation using time domain integration, Horowitz-Hubbeny 
method and exact analytical formula are included. 

then the general solution is 

Z(C) = A 1 C(1-C) /3 2^1(0,6,0,0 • (47) 

To take into account the boundary condition at spatial infinity ( = 1, we can use the 
Kummers formula to connect hypergeo metric functions at ( = with that at £ = 1, 

r(c)T(c — a — b) 

2 F 1 (a,b,c;() = ^(0,6,0 + 6-0,1-0 

r(c — a)T(c — b) 

+ (1 - O ^^lll^^ F^ ~ c - 6, c - a - 6 + 1, 1 - z) .(48) 
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2 4 6 8 10 

d 



Figure 6: Dependence of scalar quasinormal frequencies with the dimensionality of Lifshitz black 
holes. The parameters of the solution are M = m = L = k = l. 




2 4 6 8 10 

k 



Figure 7: Dependence of scalar quasinormal frequencies with the wavevector k for Lifshitz black 
holes of various dimensions. The parameters of the solution are M = m = L = 1 . 

Using (I40l-(l45l) the above expression can be written as 

Z(0 = A ie "^(l - Cf if^'! 1 ~ b,a + b-c,l-() 

L (c — a)l (c — b) 

+ A ie ^(l - 0^ T{C) 5^t' C \ Hc - a, c - b, c - a - b + 1, 1 - C) .(49) 
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Figure 8: Dependence of scalar quasinormal frequencies with the field mass m for Lifshitz black 
holes of various dimensions. The parameters of the solution are M = k = L = 1 . 




2 4 6 8 10 



Figure 9: Dependence of scalar quasinormal frequencies with the black hole mass M for Lifshitz 
black holes of various dimensions. The parameters of the solution are k = m = L = 1. 

Now, at ( = 1 we have, 

~, . -n W1 ^ r(c)r(c-a-6) p_ p T(c)T(a + b-c) 

Z(z 1) = Ax(l - Cr r(c _ fl)r(c _ 6) + Ai(l - C) » r(fl)r(6) ■ (50) 

The first term in the above expression tends to zero at spatial infinity, thus is consistent 
with the Dirichlet boundary condition in this limit. In the following, we exclude terms that 
are present in hypergeometric functions considering only non integer values of c — a — b. The 
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power of 1 — C in the second term diverges for real scalar field masses and we need to impose 
that the quotient of the Gamma functions in this term vanishes, allowing Dirichlet boundary 
conditions to be satisfied. This is true only if either T(a) or T(b) diverges, what leaves us 
with two conditions that can be used in order to determine the quasinormal frequencies, 
that is, either a = —n or b = —n, where n = 0, 1, 2, .... Using (1391) . (14"U|) and (14"51) we find for 
the quasinormal frequencies the exact result 

i {L 2 k 2 + M [(2n + l)VD 2 + Am 2 L 2 + An{n + 1) + L 2 m 2 + 5] } 
Un ~'~L [2(2n + l) + VD 2 + 4m 2 L 2 ] 

From the condition b = —n we obtain the same result for the quasinormal frequencies. The 
above expression indicates that for real scalar field masses all quasinormal frequencies are 
pure imaginary and negative, for all dimensions. Thus, higher dimensional Lifshitz Black 
Holes are stable under real mass scalar fluctuations. 

For imaginary masses, we can distinguish two cases. Putting m 2 = — |m 2 |, if (y) — \m 2 \ > 
0, we have f3 > and y — (3 > 0. Thus, the powers of 1 — C in the two terms of equation 
(15 Op vanish at spatial infinity for all values of a, b and c, as required by Dirichlet boundary 
conditions, leading to a continuous spectrum of quasinormal frequencies. For imaginary 
masses (~j — \m 2 \ < and there are waves present at the spatial infinity, that frustrate 
the fulfillment of the quasinormal boundary conditions, thus we exclude this case from our 
analysis. 

Using (I5ip it is easy to confirm the results of the last section, as is seen in Figures [6] to 
[9J in which we show the dependence of scalar quasinormal frequencies with the parameters 
involved in the perturbation problem, D, k, m and M. Again, the propagation of scalar 
waves outside higher dimensional Lifshitz black holes is a very simple phenomenon. 

VII. CONCLUDING REMARKS 

We have studied the propagation of a minimally coupled scalar field outside higher di- 
mensional asymptotically Lifshitz black holes, obtained from a theory that includes the more 
general quadratic corrections to Einstein gravity. We have solved the equation of motion 
analytically for all dimensions and performed its numerical study, obtaining the time domain 
evolution and the quasinormal frequencies. 

The quasinormal stage at late times is characterized by the absence of proper oscillations 
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and the field is purely damped. We have negative purely imaginary frequencies, implying 
complete stability of the black holes under scalar fluctuations. This picture is the same for 
all dimensions independently of the variation of parameters that describe the background 
and the perturbations, such as the black hole mass M, the scalar field mass m and the 
wave-vector of the propagating fluctuation. As all these parameters increase their values, 
Lifshitz black holes becomes more rigid, with a quicker decay of the probe field outside the 
event horizon. 

The analytical result for the quasinormal spectrum shows us that the fundamental mode 
for any dimension is purely imaginary. Concerning the gauge/gravity duality, such a result 
means that the holographic field theory has a small thermalization timescale. In the z = 3 
NMG case, the same result has been found |22j and we conjecture that the purely imaginary 
quasinormal frequencies are a universal propriety of the black holes with Lifshitz scaling 
symmetry. 
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